Two rods confined by positive plates: Effective forces and charge distribution profiles 
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Monte Carlo simulations were employed to study two negative rods confined between two parallel 
plates. The system is immersed in a 1-1 restricted primitive model electrolyte. Ion distributions 
and forces per unit of area (pressures) on rods and plates are accessed by sampling the NVT 
ensemble. Pressures are analyzed by means of their corresponding electric and contact (depletion) 
contributions. This was done for several charge distributions on plates, plates charge densities, and 
plate-plate surface separation distances. We found an enhancement of the inherent repulsive rod-rod 
effective force when uncharged plates are the confining species. On the contrary, a strong decrease of 
the rod-rod effective repulsion was obtained for positively charged plates. Moreover, attraction was 
also found for plates having a charge equal to that of fully charged bilipid bilayers. These results 
agree with DNA-phospholipid complexation experiments. On the other hand, for a model having 
the plates charges fixed on a grid, very long range rod-rod sin-like effective forces were obtained. As 
explained in the text, they are a consequence of the rod-plate double layer coupling. 

PACS numbers: 61.20.Qg,68.08.-p 



I. INTRODUCTION 

In soft matter science, confinement is known to induce 
ordering, produce density fluctuations, provoke dielectric 
changes on liquids, affect diffusion processes, promote 
layering, enhance miscibility in polymer blends, and yield 
anomalous attractions between like charged colloidal par- 
ticles JlH4j . This last effect was observed experimentally 
by accessing the pair correlation function by means of 
video microscopy and optical tweezers @. Additionally, 
biological setups as those formed by plate-like charged 
lipid bilayers and rod-like DNA molecules show that, un- 
der certain conditions, confining bilayers induce a dense 
parallel 2D array of like-charged DNA molecules [(j-fioj . 
Hence, it seems clear that confinement and self assem- 
bling are strongly linked. 

The extensively used theory to explain interactions 
among charged colloidal particles is the one develope d by 
Derjaguin, Landau, Verwey, and Overbeek (DLVO) [llj . 
This is based on the Poisson-Boltzmann theory which 
deals with the electrostatic interaction among charged 
colloids and ions. This theory predicts that like charged 
colloids should repel each other following a screened 
Coulo mb p otential. There is, however, experimental ev- 
idence [l2[ and theoretical prove [l3| that under certain 
conditions two like charged macroparticles attract each 
other in bulk, and so, this mean field theory cannot pre- 
dict this unintuitive behavior. On the other hand, there 
is also experimental evidence that like charged colloids, 
which do not attract in bulk, attract for the same con- 
ditions but under confinement Again, theories based 
on Poisson-Boltzmann cannot capture this peculiar be- 
havior. Additionally, theories based on integral equa- 
tions, which were success to predict the attraction of like 
charged colloids in bulk 0, [l5[ , are of hard implementa- 
tion for complicated geometries such as these produced 
by planar confinement. This is, at least partially, why 



there are to our knowledge not theoretical studies at- 
tempting to explain such anomalous behavior. 

Since geometry is not a problem for implementing sim- 
ulations, this turns a natural tool for studying a system 
counting on confining charged plates, confined macropar- 
ticles, the corresponding counterions, and added salt. 
This is why we employ them to reveal the effect that 
confinement plays on like charged colloids effective inter- 
actions. Additionally, since lipid bilayers-DNA experi- 
ments show transitions from loose to compact DNA 2D 
ensembles by modulating the bilayers charge density, we 
would like to focus on modeling this particular system. 
For that purpose, two DNA molecules are modeled as two 
confined like charged negative rods, and the lipid bilayers 
as the confining positive charged plates. Hence, different 
plate-plate distances, charge on plates, and charge distri- 
butions on plates are studied, for a fixed 1-1 electrolyte 
concentration. 

The paper is structured as follows. Section U is this 
brief introduction. Section [H] gives details on how the 
simulations were performed. Results are presented in 
section IIII1 and section IIVI summarizes the main con- 
tributions and tackles some conclusions. 



II. SIMULATIONS 

Monte Carlo (MC) simulations were performed to 
study a system counting on two infinitely large, parallel, 
negative rods confined by two parallel, positive plates of 
finite thickness. These hard macroparticles are immersed 
in a 1-1 restricted primitive model (RPM) electrolyte. 
The simulation box has sides lengths L X ~L Z — 200 A, 
and L y = 125 A. The origin of coordinates is set at the 
box center. Plates are located parallel and symmetric 
to the z=0-plane, with a surface-surface separation dis- 
tance r. Rods axes are placed parallel to the y axis at 
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x=±(R + t/2),z=0, being R the rods radii and t the 
rod-rod surface distance. 

For simplicity, plates, rods, and solution have the same 
dielectric constant. RPM consists of hard spheres with a 
centered point charge, such that their electrostatic inter- 
action is 

U E ir tj ) = ^p, (1) 

being /3=1/}cbT, Ub the Boltzmann constant, T=298K 
the absolute temperature, Zi and Zj the valences of sites i 

and j, £ B =^r=7.14:A the Bjerrum length, e=78.5 the di- 
electric constant, and the interparticle distance. The 
electrostatic interaction among macroparticles and be- 
tween macroparticles and ions are also given by equation 
[T]with i and/or j running over the macroparticles sites. 

We assigned to plates a charge density of 
(To=0.229C/m 2 (+e per 70 A ) for each surface, 
in correspondence with a fully charged cationic lipid 
membrane [l6j . To study the effect of varying the plates 
charge density, for some systems we assigned a fraction 
of this charge, expressed as a percentage. Three ways 
of distributing the charge on plates were explored. One 
consist of placing +e-charges in the center of 5 A hard 
spheres, which are confined by the plates boundaries and 
are moved by following the MC criteria. The hardcore 
interaction of these 5 A spheres is switched on just 
among themselves. The second way is by placing the 
+e-charges over triangular grids, which are randomly 
moved in the x and y directions to mimic a continuous 
charge distribution. These grids are placed on the 
surfaces of the plates. The third way is by placing the 
+e-charges over triangular grids, as in the previous case, 
but fixing them in these positions. On the other hand, 
rods have an — e-sitc each 1.7 A, as DNA molecules [l6j . 

Rods have a radius of i?=10.5 A, consistent with a 
hydrated DNA molecule. Plates are 5 A thick to allow 
fluid-fluid correlations between the interplate and bulk 
sides of the plates These correlations are expected 
in real 40 A thick bilayers, due to their low dielectric con- 
stant. We fixed the concentration of the 1-1 electrolyte to 
p s =0.1 M, and assigned its diameter to a=4.25 A, unless 
otherwise indicated. 

Additional anions are added to make the system elec- 
troneutral. Periodical boundary conditions are set for 
all directions. The Ewald summation formalism is em- 
ployed to deal with Coulomb interactions [l?], EH. To 
access all phase space volume, i. e., to allow for ion in- 
terchange between the confined and unconfined regions, 
movements having a large maximum displacement are 
done. Figure [T] shows a front view snapshot of an equi- 
librated system having r — 2R=0, t—a, and a 100% of 
plates charge. Here, the plates charge was distributed by 
the first method, blue means negative and red positive, 
and sizes are scaled except the 5 A plate charges which 
are represented smaller. 

Electrostatic contributions to the forces acting on the 
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FIG. 1. Snapshot of an equilibrium configuration of a system 
counting on two rods sandwiched by two positive, confining 
plates. 



macroparticles (rods and plates) are obtained by 

Fei = <EE- vc/ £M>, ( 2 ) 

i j 

where i runs over the sites of the reference macroparticle 
and j runs over all other sites. The contact (depletion) 
contribution is obtained by integrating the ions contact 
density,p(s=cte) , 

F c = -k B T J p(s = cte)nds, (3) 

where s refers to the macroparticle surface and n is a 
unit normal vector. These two contributions to the force 
are interdependent. 



III. RESULTS 

Before presenting the results obtained for confined 
rods, it is convenient to have in mind the most important 
features of the interaction between two parallel rods im- 
mersed in a 1-1 RPM electrolyte. This is why we show 
in figure [2] the effective forces per unit of rod area (pres- 
sure) as a function of the distance between rods surfaces, 
t, for two isolated rods immersed in a 1-1 electrolyte 
with a bulk concentration of 0.1 M. Note that system 
symmetry make the rod-rod forces to have only compo- 
nents on x. These components are shown in figure [5J 
where a positive sign means repulsion. There, the total 
force results from the sum of two interdependent con- 
tributions: electric and contact (depletion). It is seen 
that both contributions are always repulsive and so it 
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FIG. 2. Forces on rods per unit of rod area (pressure) as a 
function of t for two unconfined rods. 




FIG. 3. p e i(x,z) as a function of h. Darker tones indicate 
larger absolute values of p e i(x,z), while red means positive 
and blue negative. 



FIG. 4. Forces on rods per unit of rod area (pressure) as a 
function of h when highly charged confining plates are present. 

a) and b)) figure [3] shows large p e i values at (x = 0, z = 
±yjRl - {R + t/2) 2 ), being R a =R + a/2. This indicates 
the formation of two lines of counterions along the rods 
dumbbell. These counterions pressure on both rods pro- 
ducing strong repulsive contact contributions, which are 
not compensated by the external counterions [l9| . Hence, 
as t increases, the lines of counterions approach each 
other producing larger contact contributions on the x 
axis. This explains the increasing trend of the contact 
contribution with t for t < a. For t > a, the peaks at 
(x = ±t/2,y = 0) decreases with t (figures [3] c) and d)) 
and so it does the contact contribution. This produces 
the maximum of the contact component at t = a, as 
well as the maximum of the total pressure. Finally, it 
should be mentioned that the work per unit of charged 
site necessary to move the rods from i=4a to t—0 is 2.13 
K~bT / site. 

A. Influence of confinement 



is the net force. The electric contribution is ever de- 
creasing, since likely charged rods inherent repulsion de- 
cays with distance and since more counterions interpose 
between them as t increase. When rods double layers 
turn independent the one another, both contributions be- 
come cero, as expected. The contact contribution, how- 
ever, behaves differently. It increase with distance for 
t < a, and decays for t < a, showing a maximum for 
t = a. To understand this behavior, one should pay at- 
tention on how the charge distribution profiles around 
the macroparticles, p e i(x,z), changes with t. This is de- 
fined by pei(x, z)=z + g + (x, z;r,t) + Z—g—(x,z;T,t), be- 
ing g+(x, z;r,t) (g~(x, z;r,t)) the cation (anion) distri- 
bution profiles, given r and t. This case corresponds 
to t — y oo (no plates in the system). Figure |3J shows 
p e t(x,z) for t=0, 0.75a, 1.5a, and 2.25a. For i<a (plots 



Keeping all the same conditions fixed but including 
the plates in the system, we studied how confinement 
influences the rod-rod effective interaction. As a first ex- 
periment, we considered the plates to be fully charged 
(100%), having the charge distributed according to the 
first method (charge sites are confined by plates bound- 
aries but are moved by following the MC criteria). This 
model may describe the distribution of phospholipid 
charged heads inside the bilayer (plates). 

The forces obtained by means of this experiment are 
shown in figure 01 They strongly contrast with the pre- 
viously shown in many aspects. To begin with, all con- 
tributions are very close to cero for all distances, t. In 
particular, the total force is practically negligible due to 
the counterbalance of a small repulsive electric contri- 
bution and a small attractive contact contribution. This 
leads to a work per unit of charged site necessary to move 



4 




FIG. 5. p e i{x,z) for t=0 and h= 0, 0.75a, 1.5a, and 2.25a. 
Darker tones indicate larger absolute values of p e i(x, z), while 
red means positive and blue negative. 



the rods from i=4a to t=0 of 0.01 K B T/site. In addi- 
tion, it seems to be attractive in the region 1.25a<t<2a, 
and clearly positive for t<a. In fact, we obtain a very 
small but negative work per unit of rod site necessary 
to move the rods from i=4a to i=1.25a, i. e., -0.0005 
KBT/site. Hence, strong changes in the charge distri- 
bution profiles around the macroparticles are expected. 
They are shown in figure [5] for h~0, 0.75a, 1.5a, and 
2.25a. As can be seen, the charge density surround- 
ing the rods, which used to be positive, is now nega- 
tive. It should be noted that this case corresponds to 
the one shown in figure [H where just a cation is seen 
in the confined region for this particular configuration. 
Hence, cations are expelled out from the sandwiched re- 
gion, in favor of anions, which are the plates counterions. 
This large anion concentration in the interplate region 
counterbalances the rod-rod electrostatic repulsion, pro- 
ducing a remarkable decrease of the electrostatic contri- 
bution. On the other hand, the large cation peaks at 
(x = 0, z = ±yi2^ — (R + t/2) 2 ), which were responsi- 
ble for the large contact repulsive contribution, are not 
present any more, and so it is the contact repulsive be- 
havior. Moreover, the anions peaks that grow outside the 
rods and close to the pl ates overcompens ates the anion 
peaks at (x = 0, z = ±V#a — (R + t/2) 2 ) producing an 
attractive contact contribution. These outside peaks are 
seen in figure [5] for all distances t. It is also seen that 
the inner peaks at (x = 0, z = ±y/ R 2 — (R + t/2) 2 ) are 
smaller for figures [5] a) and b). In fact, figure [5] a) shows 
very narrow peaks at (x — 0,z = ±y / R 2 — (R + 1/2) 2 ) 
indicating that anions can barely enter in these inter rod 
regions. 

Figure [5] also shows large interplate ionic depletion re- 




FIG. 7. pei{x,z) for t=a and h= (a), 0.75a (b), 1.5a (c), 
and 2.25a (d). Darker tones indicate larger absolute values of 
p e i(x, z), while red means positive and blue negative. 



gions due to the presence of the rods. As expected, these 
depletion regions unbalance the contact forces inside and 
outside the plates, producing an attractive net contact 
force per unit of plate area of -6.72 MPa. This force 
is practically independent of the rods positions (the in- 
formed value is the average) . The electric contribution to 
the force on plates is repulsive but practically cero, i. e., 
1.12 MPa. This leads to a total force per plate area of 
-5.6 MPa. 

A second computer experiment was carried out for the 
same conditions but r — 2R=^2a. Under these condi- 
tions, the forces per unit of rod area are plotted in figure 
IH1 These forces show an intermediate behavior between 
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FIG. 8. pli(x,z) for planes. Darker tones indicate larger 
absolute values of p e i(x, z). 

cases t — 2R—0 and r — 2i?— s-oo (compare figures [5] and 
[4] with figure [6]). That is, both contributions to the force 
were diminished by confinement, but in this case, the de- 
crease is not as pronounced as for r — 2R=0. Indeed, 
the work per unit of rod sites to move the rods from 4a 
to is 0.78 KbT/site, which lies in between the values 
obtained for unconfined and confined conditions. Addi- 
tionally, other features such as the contact force peak at 
t=a persist, which is a signature of the presence of rods 
counterions in between the rods, as found for r — 2i?^»oo. 
Nevertheless, the negative contact contribution for short 
rod-rod distances indicates the presence of large anions 
peaks outside the rods, as found for r — 2i?=0. 

These features are depicted in figure [7] where the charge 
density profile is shown for t=0, 0.75a, 1.5a, and 2.25a. 
There, the yellow regions close to the rods and at the mid- 
plane indicate the presence of rods counterions (cations). 
It is observed that they enter in between the rods as far 
as their surface-surface distance increase. This is similar 
to case r — 2i?— i-oo. In addition, anion peaks outside the 
rods are larger than those in between them for t=0 and 
0.75a (figures[7]a) and b)). Another important difference 
with respect to case r — 2R=0 is that interplate ionic de- 
pletion regions are smaller. This impacts on the contact 
plate-plate pressure, which in this case is -3.6 MPa. The 
electric contribution turns attractive, i. e., is -1.0 MPa, 
and the total pressure is -4.6 MPa, for t = 0. In this 
case, the attractive pressure on plates decreases with t 
reaching values close to -3 MPa. 



B. Influence of charge distribution on plates 

It was shown in the previous section that confinement 
strongly affects the interaction between two macroparti- 
cles (rods, in our case). In the present section, we would 
like to elucidate the influence of the charge distribution 
on plates on the rod-rod effective force. To that purpose, 
we should first pay attention on the positive charge dis- 
tribution on plates, p*i(x, z), obtained for case r — 2R=0, 
t = 0. Here, p* el {x,z) is defined by z*+g*+(x, z; r, t), be- 
ing z* + the valence and g*+(x, z; r, t) the distribution 
profiles of charged sites inside the plates, given r and 
t. For these conditions, p* el {x,z) is shown in figure [S] for 
the lowermost plate, where the rods are schematized for 
clarity. 

It can be seen that, driven by energy and entropy, the 
plates charges move toward the plates surfaces forming 
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FIG. 9. Forces on rods per unit of rod area (pressure) as a 
function of h for a continuously like surface charge distribu- 
tion on plates. The lines corresponds to the results presented 
in figure [4] 



almost homogeneously distributed charge surface densi- 
ties, which are not very influenced by the rods positions. 
This is maybe the most important conclusion, since it 
makes us think that a homogeneous surface charge distri- 
bution may also model this system. On the other hand, 
it is also true that a small deviation from an homoge- 
neously surface charge distribution is obtained. That is, 
charge peaks close to the rods are observed, which in turn 
produce relative depletion regions at their sides. These 
features are highlighted in figure [8] 

Hence, let us test a simpler model where the posi- 
tive charge on plates is distributed on triangular grids 
at both plates interfaces, which are frequently, indepen- 
dently, and randomly moved in the x and y directions, 
to mimic homogeneous surface charge distributions. The 
forces on rods obtained by using this model are presented 
in figure [9] and compared with those obtained by mov- 
ing the charged plate sites by following the MC criteria. 
As can be seen, in general, a very similar behavior is 
obtained. That is, both entropic contributions are neg- 
ative, monotonously increasing functions of t, yielding 
cero for t>2.5a, and both electric contributions are posi- 
tive, ever decreasing functions of t, producing cero values 
for t>2.5a. Additionally, both total contributions yields 
practically cero values for all t, being positive for t<a and 
negative for 1.25a<i<2.25a. The main difference seems 
to be the smaller value of the total force at rod-rod con- 
tact obtained for the continuously like surface charge dis- 
tribution model. On the other hand, the obtained value 
of the total plate-plate force is -5.5 MPa, which com- 
pares well with that obtained for movable plates charged 
sites (-5.6 MPa). Hence, we may conclude that, at least 
for these conditions, the discrete nature of the bilay- 
ers charge should not strongly affect the forces on DNA 
molecules. 

Leaving aside the small differences between the previ- 
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FIG. 10. Forces on rods per unit of rod area (pressure) as a 
function of h for a continuously like surface charge distribu- 
tion on plates. 



ous studied models, we still would like to know whether a 
structured fixed grid of discrete surface charge affects the 
rod-rod effective interaction. This, of course, is not ex- 
pected to be a good model for a lipid bilayer, but may be 
a model of certain type of nanostructured layers. Results 
are shown in figure fTTJl For very small t, they are similar 
to those shown by figure^ i. e., there is an attractive 
decaying contact contribution and a repulsive decaying 
electrostatic contribution. For t> 0.5a, however, the con- 
tact contribution crosses over yielding positive values. In 
addition, it produces a maximum at £=1.25a. Moreover, 
it produces positive maxima and negative minima at ev- 
ery 7.8 A, showing a sin-like behavior. Note that these 
maxima and minima reach similar absolute values, inde- 
pendently of how large the rod-rod distance is. In other 
words, this behavior is very long ranged and so, it seems 
to be a consequence of the rod-plates interaction, indeed. 
On the other hand, the electric contribution shows also 
a sin-like behavior, once the rods are far enough the one 
another. Its amplitude is much smaller than the contact 
contribution, but seems to be in phase with it. 

The rod-rod effective forces must be related to the 
charge distribution profiles of the interplate region. This 
is shown in figure [TT] for h=0, 0.75a, 1.5a, and 3a. In all 
plots of figure [TT] are seen peaks and valleys of the charge 
distribution close to plates. Naturally, peaks coincide 
with the locations of the positive charge of plates and 
valleys are between peaks. Hence, the peak to peak dis- 
tance is equal to the x distance between adjacent charged 
plates sites, A=7.8 A. Note that the locations of the in- 
ternal charged sites of both plates are the same by con- 
struction, and so, internal peaks and valleys of the charge 
distribution profile of one plate face those of the other. 
This produces a reinforcement of the effects that peaks 
and valleys produce on the rod effective forces. For in- 
stance, in figure [TT] a) it is observed that the anionic dis- 
tribution at contact with the internal (closer to the other 
rod) surface of rods corresponds to valleys and outside to 




FIG. 11. p e i(x, z) for a discrete surface charge distribution on 
plates, r — 2i?=0, and h= (a), 0.75a (b), 1.5a (c), and 3a 
(d). Darker tones indicate larger absolute values of p e i(x,z). 
Red dots in plates indicate where the discrete positive charge 
is located. 



peaks, explaining the attractive contact force at t=0. In 
figure [TTIb). the inside and outside charge distribution at 
rod surfaces do not correspond to peaks or valleys, and 
hence, contact forces correspond neither to maximum nor 
to a minimum (i=0.75a). For <=1.5a a maximum is ob- 
served for the contact force, which is explained by the 
ionic peaks at contact with the internal surface of rods 
(see figure [TTIc)). and the valleys at contact with the ex- 
ternal rods surfaces. The opposite situation is found for 
£=3a, where peaks are found in contact with the outside 
of rods and valleys at the inside (see figure[TT]d)). Hence, 
it turns clear that this apparent long range rod-rod in- 
teraction is nothing but a short range rod-plate double 
layer coupling effect. 

It must be noted that this kind of effects are a function 
of \f{R 2 — {R— a/2) 2 )/A and of the relative position be- 
tween the charged sites of one plate and the other. For in- 
stance, if the charged sites of one plate and the other had 
not been face-to- face but unphased a distance A/2, none 
of these effects would have occurred. None of these ef- 
fects would have happened for y/{R 2 — (R— a/2) 2 )/A=0.5 
either. On the contrary, these effects are expected to be 
strong for \/{R 2 — (R — a/2) 2 )/A=l. Our system has 
^{R 2 - (R- a/2) 2 )/A= 0.82 and has the charged sites 
of one plate faced to the charged sites of the other, and 
this is why the effects are clearly seen. Another thing 
to note is that the x position of one plate with respect 
to the other controls the phase of the two lines of peaks 
and valleys of the charge distribution profiles close to 
plates, which in turn controls the apparent rod-rod long 
range interaction. Hence, by changing the plates relative 
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FIG. 12. Forces on rods per unit of rod area (pressure) as a 
function of h for plates having a 50 % charge density. 



positions the rod-rod separation distance should change 
accordingly. This property may be of practical interest 
for nanostructure manipulation. 




FIG. 13. p e i(x,z) for 50 % charged plates, T — 2R—0, and h= 
(a), 0.75a (b), 1.5a (c), and 2.25a (d). Darker tones indicate 
larger absolute values of p e i(x, z). 



C. Influence of plates charge density 

DNA-lipid bilayer lamellar phase ensembles show a 
transition from condensed to loose DNA 2D arrays mod- 
ulated by the rate of cationic (di-oleoly thrmethylam- 
monium) to neutral (dioleoyl phosphatidylcholine) lipid 
molecules, which form part of the confining bilayers. Ex- 
periments show that for a fully charged membrane, i. e., 
no neutral lipid added, DNA form a condensed arrenge- 
ment which has a DNA-DNA distance of 24.5 A. This 
distance enlarges to 57.1 A for 57% of the full charge 
density. Hence, we find interesting to study the rod-rod 
effective interaction for smaller plates charge densities. 

Results for half the full charge density of plates are 
shown in figure [12] The total effective rod-rod force is 
a positive monotonically decreasing function of the rod 
surface distance, i. e., its always repulsive. Hence, the 
work per unit of rod site to bring them in touch is 1.15 
KgT '/ site. It should be noted that this value is relatively 
similar to that obtained for fully charged membranes but 
for a larger interplate distance (0.78 KbT/ site). In fact, 
all force components of figure [12] are very similar to those 
shown in figure [6] On the other hand, the pressure on 
plates is -5.1 MPa for the rods at contact, and decreases 
its absolute value to -4.5 MPa for larger rod-rod dis- 
tances. In this case, the electric component is the larger 
contribution being -4.8 at i=0, and decreasing to yield 
-3.2 for t=4a. 

The corresponding charge densities are shown in figure 
[131 As expected, much lower values of the charge den- 
sity at contact with the plates surfaces are obtained, as 
a consequence of the lower charge of the plates. Hence, 
the interplate region becomes less dense in anions (plates 
counterions) and positive charge densities are observed 
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FIG. 14. Forces on rods per unit of rod area (pressure) as a 
function of h for uncharged plates. 



close to the rods surfaces. In particular, cations are 
allowed to fill the inter rod regions at (a; = 0, z = 
±yi£| — (R + t/2) 2 ), which are responsible for the con- 
tact characteristic peak at t = a (similar to figures [5] and 
[5]). Another consequence of a less crowded anionic in- 
terplate is the increase of the rod-rod effective repulsion. 
That is, the less counterions in the interplate, the less 
compensated the inherent rod-rod repulsion becomes. Fi- 
nally, other point to note is that the charge densities at 
contact with the outer plates surfaces becomes clearly 
dependent on the x coordinate. Since in this case the 
charge on plates do not overcompensate that on rods, 
the presence of the rods strongly affects the outer charge 
distribution. 

When confinement is produced by uncharged plates 
the rod-rod total effective force becomes very repulsive, 
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FIG. 15. p e t(x,z) for uncharged plates, r — 2R—0, and h= 
(a), 0.75a (b), 1.5a (c), and 2.25a (d). Darker tones indicate 
larger absolute values of p e i(x, z). 

mainly due to the increase of the contact contribution. 
This is clearly seen in figure [T4J where the large peak of 
the contact contribution at t = a makes the total pressure 
to reach values close to 20 MPa. The general trend of the 
contributions is, however, similar to that already shown 
for unconfined conditions. In addition, the obtained pres- 
sure for plates is 2.4 MPa for t = 0, which decrease for 
increasing t yielding 0.4 MPa for t = 4a. Obviously, the 
total pressure coincides with the contact contribution, 
since a cero electric contribution is obtained. Due to the 
fact of obtaining large positive plate-plate pressures, a 
configuration where rods are squeezed by plates is, in 
this case, not expected. 

Confinement by uncharged plates also affects the 
charge distribution profiles. The fact that no ions can 
attach to part of the rods surfaces due to the plates ex- 
cluded volume, makes the system react by pushing more 
ions to the available surfaces. In particular, those pore- 
like regions (large surface/volume rate), as those located 
at (x = Q,z = ± ^JR\ - {R + t/2) 2 ), turn crowded by 
cations (rods counterions). Again, these cations exert 
large repulsive contact contribution on rods, explaining 
the huge effective force peak at t = a. 



IV. CONCLUSIONS 

Monte Carlo simulations were performed to study a 
system counting on two charged, parallel rods, and two 
confining charged, parallel plates. The system is im- 
mersed in a 1-1 restricted primitive model electrolyte. 
Effective forces acting on the macroparticles and the cor- 
responding charge density profiles were accessed for dif- 
ferent confinement conditions, i. e., the influence of the 
amount of charge of the confining plates, its distribution, 



and the plate- rod distance were studied). Both situa- 
tions, an enhancement of the inherent repulsive rod-rod 
effective force, and a strong decrease were obtained for 
uncharged and charged plates, respectively. Some cases 
even show attraction, as found for highly charged confin- 
ing plates. Hence, highly charged plates promote rod-rod 
crowding, and lower charged plates induce larger rod-rod 
separation distances. This well agrees with DNA-bilipid 
experiments. We found that models having movable 
charged sites on plates and homogeneously-like charged 
surfaces lead to very similar results. 

On the other hand, a model having discrete, fixed site 
charges on plates, located forming a regular grid, induce 
a very long range rod-rod effective interaction. This long 
range interaction shows a sine-like behavior, having sev- 
eral maxima and minima regularly separated. As pointed 
out in the paper, it is a consequence of the double layer 
coupling between rods and plates, and has little to do 
with the rod-rod direct correlation. This is why we pre- 
fer to qualify it as apparent long range rod-rod interac- 
tion. Since this interaction rules the rod-rod separation 
distance and since the interaction depends on the rela- 
tive position between plates, controlling this relative po- 
sition, in turn, one controls the degree of compactness of 
the parallel stack of rods. This may be of practical use 
for nanoconstruction purposes. 

In brief, confinement and self assembling are strongly 
connected, since attraction between rods and plates de- 
pends on a delicate balance between charge and excluded 
volume of all particles in the system. Additionally, self 
assembling is linked to several interesting phenomena 
such as double layer coupling, counterion release, local 
violations of electroneutrality, and to the appearance of 
ionic depletion regions. Our belief is that neither of 
them by itself rules self-assembling, but self-assembling 
involves all of them acting interdependently, as a whole. 



Appendix A: Electric field 

We find instructive to show how the electric field looks 
like for some of the cases here studied. We restrict to 
show the t=0 configurations, which correspond to the a) 
plots of figures El El [Zl El [131 HS1 i- e., to unconfined 
conditions, confinement by 100% charged plates at con- 
tact with rods, by fully charged plates at r — 2R = 2a, 
by plates with fixed charges, by 50 % charged plates, 
and by uncharged plates, all these three last cases for 
t — 2R = 0. The electric field at a given x y position 
was calculated by simply averaging the electric force per 
unit of charge acting on charged sites located at the cor- 
responding position. This average is performed for all 
configurations. For obtaining the electric field for those 
x and y positions where charged sites cannot access, an 
extra ghost charged site is included in the calculation. 
Naturally, this extra ghost charged site is not accounted 
during the simulation. 

The obtained electric fields are shown in figurcQl)! Fig- 
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FIG. 16. Electric field for the different cases presented along the paper. Arrows indicate the field direction and darker tones 
points to larger absolute values of it. Tone changes are proportional to log(|E|) 
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ure[16ja) corresponds to the unconfined case. There it is 
seen that the electric field converges to both rods axes, 
as expected. The intensity of the electric field decays 
with the distance to the rods axes, for the region outside 
the rods, and in a more complicated way for the region 
in between them. In particular, there is a point of cero 
electric field at x — z = 0, which is noted as lighter 
tones of the surrounding arrows. 

The electric field is very disturbed by the presence of 
the confining plates, whether they are charged or not. 
This is in general seen in figures [TBI bVfh Case b), fully 
charged plates at contact with rods, shows that inside 
plates the electric field is very week, due to the fact that 
plate's charge mostly locates at the interface, as com- 
mented in section Hill This fact makes difficult to obtain 
accurate values for this region, and this is why noise ap- 
pears in the arrows directions. It should be noted that 
small positive (repulsive) values for the electric compo- 
nent of the force acting on plates were obtained under 
these conditions. This agrees with the fact of having a 
very small electric field inside the plates. In addition, and 
as expected, the electric field diverges from the plates. 
This, and the fact of having convergence at the rods axes, 
produces a very intense local electric field in between rods 
and plates. There also appears a region of cero electric 
field at z = 0, outside the rods. A final thing to note is 
that the presence of the rods clearly affects the electric 
field outside the interplate region. 

When plates are a little bit separated from the rods, 
i. e. for r — 2R = 2a (case c)), it is seen that the elec- 
tric field inside plates increases. This may suggest that 
larger positive electric force acting on plates should be 
obtained for this case. However, this force is relatively 
large but negative (attractive), as mentioned in section 
IIIII A closer inspection of figure [T5] c) reveals that arrows 
diverging from plates toward the midplane are darker 



than those diverging outside plates. This means that 
charged sites on the internal plates surfaces suffer larger 
forces than those at the outer surfaces, which in turn 
explains the attractive electric force on plates for this 
case. Additionally, it is seen that these darker arrows are 
located close to the rods, and lighter arrows are seen di- 
verging from plates for larger distances to the rods. This 
means that the electric component to the plates pres- 
sure is concentrated close to rods, suggesting that plates 
are supporting shear stress. Hence, non rigid plates 
such as bilipid membranes may deform to surround rods, 
as shown experimentally for certain ADN-phospholipid 
complexes. 

Figure [Tol d) corresponds to the case in which we fixed 
the charged sites at plates surfaces. These charged sites 
are clearly distinguished in figure |T5] d) since they act as 
divergent sources of electric field. This, in turn, produces 
the interplate charge density peaks and valeys already 
shown in figure 1111 

Finally, figures [T51 e) and f) show the results for plates 
charged with 50 % of the full charge, and uncharged 
plates, respectively. Figure e) shows that the electric 
field inside the plates points toward the inter rod re- 
gion, and that arrows that diverge from plates toward 
the plates midplane are darker than those pointing in 
the other direction. This fact also well agrees with the 
fact of obtaining an attractive electric contribution for 
the pressure acting on plates, as mentioned in section. 
Again, this case suggests that shear stress is being sup- 
ported by plates, since the electric field focus just in the 
region close to rods. On the other hand, figure [15] f) 
shows that the electric field diverges from some space 
regions which concentrate positive ions. For instance, di- 
vergences are seen at (x = 0, z = iyiJjj — (R + £/2) 2 ), 
at (x = ±25 A, 2 = 0), and at z = ±22 A. These regions 
show positive charge distributions, as clearly shown by 
figure [151 
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